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Recent progress in calculations of QED radiative corrections to parity nonconservation in atoms 
is reviewed. The QED vacuum polarization, the self-energy corrections and the vertex corrections 
are shown to be described very reliably by different methods used by different groups. All new 
calculations have recently converged to very close final values. Each separate radiative correction 
is very large, above 1 % for heavy atoms, but having different signs they partly compensate each 
other. Our results for the radiative corrections for all atoms are presented. The corrections are 
-0.54% for 133 Cs, and -0.70% for 205 T1, 208 Pb, and 209 Bi. The result for 133 Cs reconciles the most 
accurate atomic experimental data for the f 
the standard model. 

PACS numbers: 32.80.Ys, 11.30.Er, 31.30.Jv 
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I. INTRODUCTION 

It has become clear quite recently that the QED cor- 
rections give a large contribution, of the order of ~ 1%, 
to the parity nonconservation (PNC) amplitude in heavy 
atoms, which makes them important for an analysis of 
modern experimental data. The present paper summa- 
rizes recent progress in calculations of the QED radiative 
corrections to PNC. Over the last few years the experi- 
mental data of the Boulder group for the 6s — 7s PNC 
amplitude in 133 Cs, which remains the most accurate in 
the area of atomic PNC, were presumed to be able to 
challenge the standard model. The radiative corrections 
eliminate this possibility, their accurate determination 
brings the atomic experimental data for the mentioned 
transition in agreement with the standard model. The 
paper presents also a brief description of the methods 
used for calculations of the QED corrections. Some of 
them, that have been recently designed for applications 
to the PNC problem, may be of interest in other areas 
where QED corrections are studied. 

Several atoms and several transitions have been stud- 
ied experimentally in relation to PNC. As was mentioned, 
the most accurate experimental data was obtained for 
the PNC 6s - 7s PNC transition in 133 Cs. Studies of this 
transition initiated by Bouchiat and Bouchiat ll|, were 
continued and developed in Paris 0, 0, 0, IE IE 01 an d in 
Boulder @, EE El The latest work of Wood et al [H 
reduced the experimental error to 0.35% that remains the 
best accuracy in atomic PNC experiments. These results 
inspired interest in very accurate atomic PNC calcula- 
tions that, eventually, led to recognition of the fact that 
the QED corrections contribute at the necessary level 
of accuracy. Other experiments on PNC have been car- 
ried out for several transitions in several atoms. This 
includes the 6p 1/2 - 7p 1/2 transition in 205 T1 |T1IT3.IT^. 



and the 6p 1/2 - Gp 3 / 2 transition in 205 T1 111 HI Il7j. 
the 3 P - 3 Pi transition in 208 Pb EE EE El and two 
transitions in 209 Bi, 4 S^f 2 — 2 -D5/2 measured in pioneer- 
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ing works of Barkov and Zolotorev |2lL 
in El IM |22l as well as the 4 S 3/2 - 2 D, 
pi l29t l3C| . The comprehensive theoretical background 
on PNC as well as a detailed account of events which led 
to a discovery of PNC in atoms (that was also the first 
observation of PNC in neutral currents) can be found in 
the book of Khriplovich [3l|] . 

Errors of final quantities that can be extracted from 
experimental results on PNC in atoms crucially depend 
on the accuracy of atomic calculations. The most diffi- 
cult and cumbersome part of the theoretical research is 
the atomic structure calculations that take into account 
the many-electron nature of the atomic wave function 
in heavy atoms. A foundation for the modern theoreti- 
cal description of many-electron correlations in the PNC 
problem was laid out in the works [3E |3j| which also 
give references to previous papers. The calculations of 
the many-electron correlations in the 6s — 7s transition 
in 133 Cs in these works gave very close results, with the 
error in both sets of calculations estimated as ~ 1%. The 
problem was revisited in several recent publications that 
supported the results of I32L 33] improving the error. Im- 
provement of the theoretical error was initiated in Refs. 
[34| that compared a number of experimentally measured 
quantities, such as polarizabilities, dipole amplitudes and 
hyperfine constants, with theoretical predictions of Refs. 
[32I l33| . concluding that the theoretical error of these 
works was probably better than it had been anticipated, 
at the level of 0.4%. This conclusion allowed the authors 
of |2E| to reveal a deviation that existed between the ex- 
perimental data of EH and the standard model (~ 2a). 
The following theoretical works [IE HE 111 demonstrated 
that indeed, the theoretical error that originates from 
the manyelectron correlations can be reduced down to 
~ 0.5% for the 6s — 7s transition in 133 Cs, thus matching 
the experimental accuracy of Ref. El- The net result 
of the experimental and theoretical progress at this pe- 
riod of time was a pronounced deviation, of the order 
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of 2a, between predictions of the standard model and 
experimental data that was proved to be not related to 
uncertainties imposed by the complicated manyelectron 
nature of the problem. 

The progress outlined above indicated that either the 
standard model needed improvement, or some unrecog- 
nized, underestimated factors should come into play. A 
careful analysis revealed that there was an overlooked 
factor well within the framework of the standard model, 
namely the conventional QED corrections to the PNC 
amplitude. Obviously they have always been kept in 
mind, but for a long period of time they were estimated 
as negligible in the PNC problem. From the perspec- 
tive of the present day situation the fact that they play 
an important role can be considered as a general state- 
ment. Historically, however, its realization stemmed from 
several works devoted to different, separate phenomena 
that contribute to the QED corrections. Derevianko [3^ 
found that the Breit corrections to electron-electron in- 
teraction in heavy atoms are much larger than it had 
been anticipated, giving —0.6% for the 6s — 7s PNC am- 
plitude in 133 Cs. This result was obtained numerically, 
and later confirmed in a number of papers [sH . l39i Efij 
both numerically and analytically. The account of the 
Breit corrections put the experimental data of much 
closer to the standard model, but this was only the begin- 
ning of an intriguing zigzag road of research. An analysis 
of the Breit corrections by Sushkov |4£j revealed that one 
can expect the QED vacuum polarization to be as impor- 
tant as the Breit corrections. Johnson et al |4j for the 
first time proved that indeed, the vacuum polarization 
is very important. For the 6s — 7s PNC amplitude in 
133 Cs atom they found that it gives 0.4%, as confirmed 
in |36l l42l |43| . From a purely pragmatic, numerical point 
of view the vacuum polarization mostly compensated for 
the contribution of the Breit corrections, bringing the ex- 
perimental data for the 133 Cs atom back in contradiction 
with the standard model. 

The most difficult and, correspondingly, challenging 
part of the QED corrections to the PNC amplitude is 
presented by the QED self-energy corrections. We are 
using this term here in a sense that embraces the ver- 
tex corrections as well. One could have anticipated these 
corrections to be large from the very beginning. The self- 
energy corrections are known to give a large contribution 
to a number of phenomena in atomic physics, including 
the Lamb shift, the hyperfine interaction, and the en- 
ergy shift due to the finite nuclear size (FNS). The latter 
phenomenon will play a particular role in the following 
discussion. There is no physical reason indicating that 
the PNC amplitude should be any different. The sign 
of the self-energy corrections is also predetermined. The 
QED self-energy is known to be mostly repulsive (this is 
definitely true for the S1/2 state that should play a very 
important role) in contrast to the attractive vacuum po- 
larization. Therefore the self-energy corrections should 
make the electron wave function smaller at the nucleus 
and, consequently give a negative contribution to the 



PNC amplitude. 

However, the aforementioned reasons were confronted 
by counter-arguments indicating that the self-energy cor- 
rections should be small. One reason, briefly mentioned 
in |4l|. relies on the fact that in the plane- wave approxi- 
mation the self-energy corrections are small, of the order 
of 0.1% for 133 Cs atom, as was found by Marciano and 
Sirlin |44| and Lynn and Sandars [45| . One can argue that 
account of the nuclear Coulomb field should not produce 
any dramatic enhancement, and the self-energy correc- 
tions should, probably, remain small. The argument has 
its merit, it is known to be correct for the case of the vac- 
uum polarization, where higher-order corrections are sup- 
pressed, see section [H] However, the self-energy correc- 
tions behave differently. The well-studied case of the hy- 
perfine interaction is an example. It is known from both 
analytical and numerical calculations, see Refs. 0, 
and references therein, that in this problem the manifes- 
tations of the Coulomb field in the self-energy corrections 
in heavy atoms are more pronounced than a contribution 
of the plane-wave approximation, the latter one is given 
by the known Z-independent Schwinger term q/27r, while 
the former rise quickly with the nuclear charge |76j . 

Another case, that looked even stronger, against a pos- 
sible role of the self-energy corrections in the PNC prob- 
lem was made by Milstein and Sushkov in Ref . ■ They 
presented results of direct analytical calculations claim- 
ing that the self-energy corrections are positive and small, 
concluding that the experimental data on the 6s— 7s PNC 
amplitude in the 133 Cs atom are in contradiction with the 
standard model. The impact of this work was strong, for 
a long period of time it was considered as a reliable argu- 
ment indicating that the experimental data can challenge 
the standard model. The opposite opinion expressed in 
Ref. j3^|, which presented a model for rough estimation 
of the self-energy corrections, was not considered as con- 
vincing, probably because the proposed model neglected 
the vertex corrections and was not gauge invariant. 

The analytical calculations of Ref.^2] seemed to prove 
that the self-energy corrections are small and positive, 
while some physical arguments mentioned above indi- 
cated otherwise. To find the way out of this controversy 
we proposed in Ref. 0] an approach that avoided dif- 
ficult analytical calculations. This work derived a new 
equality, we will call it the chiral invariance identity due 
to reasons explained in section IIII Al It expressed the 
self-energy corrections for the PNC amplitude in terms 
of the self-energy corrections to the energy shifts due to 
FNS. The latter ones were calculated numerically in Refs. 
|49t l50L lo li 152 . l53l l54| . Thus the found equality allowed 
one for the first time to find the accurate self-energy cor- 
rections for the PNC amplitude for heavy atoms. As 
one could have anticipated, they were large and nega- 
tive. In particular, for the 6s — 7s transition inl33Cs we 
found —0.73(20)% that brought the experimental data of 
Wood et al in agreement with the standard model. The 
uncertainty of this result was mainly due to errors of the 
numerical results of Ref. 52] that we used as input data. 
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The most pronounced among them is the error for the 
self-energy corrections to the energy shift induced by the 
FNS in the p\/2 wave. The authors of (52], as well as 
others in the relevant papers above, estimated their re- 
sults as very reliable for heavy atoms Z ~ 90, and less 
accurate for lighter atoms, making numerical calculations 
below Z — 55 difficult. 

In order to verify and refine our results, there were per- 
formed analytical calculations of the leading, linear term 
in the Za-expansion for the self-energy (and vertex) cor- 
rections for the PNC amplitude [55l | . The result of this 
work confirmed that of |4g, the self-energy corrections 
were found to be large and negative, in clear contradic- 
tion with Ref . 42] . Combining the data of , which is 
accurate for heavy atoms, with the analytical data of [55| 
that is reliable for light atoms, Ref. [55| found the self- 
energy corrections for the PNC amplitude for all atoms. 
This analysis also allows one to reduce the impact of the 
above mentioned errors of the data extracted from [H^ . 
The result for the 6s— 7s transition in 133 Cs was —0.9%, 
in good agreement with —0.73(20)% of |^. These two 
works, using different techniques, unambiguously demon- 
strated that the experimental data of Wood et al 
support the standard model. 

The contradiction that existed between this statement 
and conclusions of Ref. 0] wa s, fortunately, eradicated 
by the the work of Miistcin et al 5j| that appeared 
shortly after 1551 . Milstein et al arrived at results that are 
close to ours |48ll55j. For the 6s — 7s transition in 133 Cs 
they found the self-energy corrections to be —0.85%. A 
careful comparison of the results of the two groups pre- 
sented in Ref. |57j | revealed a similar agreement for all 
values of Z. In this latter work it was also proposed that 
the remaining small discrepancy is due to limitations of 
accuracy in calculations of Milstein et al that are based 
on Za-cxpansion (Za = 0.4 for Cs). This conclusion 
has recently been proved true by Milstein et al who re- 
fined their calculations in Now, for a vast area of 
the nuclear charges Z the results of Milstein et al for the 
self-energy corrections agree with ours very closely, the 
remaining deviation is below 9%. The deviation mani- 
fests itself more prominently for large Z , and we believe 
that it is still due to those terms of the Za-expansion 
that remain unaccounted for by the analytical calcula- 
tions of [Hsj . Thus, using several different techniques, 
the two teams arrived to one and the same result for 
the self-energy corrections. The fact that the results of 
Milstein et al converged to ours after a period of initial 
deviation makes the final results even more convincing. 

The plan of this paper is to discuss in some detail the 
radiative corrections, summarizing all available methods 
and data. As an introduction to this study let us first 
present very briefly several known facts about relativistic 
corrections in general. The main relativistic correction is, 
obviously, the correction that is due to the relativistic na- 
ture of the Dirac equation that governs the electron prop- 
agation. It is usually taken care of in the initial approx- 
imation that relies on the Dirac-Hartree-Fock approach, 



using then the machinery of the many-body theory to 
account for the electron correlations. The next necessary 
step is to include the relativistic correction to the photon 
propagator that describes the electron-electron interac- 
tion. In the nonrclativistic approach this propagator is 
approximated by its Coulomb part. Relativistic effects 
can be taken care of perturbatively, using the Breit cor- 
rections for this propagator. This approximate approach 
is convenient because the Breit correction can also be 
incorporated in the initial Hartree-Fock mean field and 
then dealt with by conventional many-body methods. At 
the same time this approximation provides sufficient ac- 
curacy. 

The next layer of corrections originates from the QED 
radiative corrections. Numerically, as well as parametri- 
cally, they may be comparable with a contribution of the 
Breit corrections, but it is convenient to consider them 
separately. The radiative corrections in the lowest order 
of the perturbation theory are divided into two classes. 
One of them constitutes the QED vacuum polarization, 
another one consists of the self-energy and the vertex cor- 
rections. These latter two corrections are often grouped 
together and called the e-line corrections. 

The discussion below focuses entirely on the QED ra- 
diative corrections. This presumes that every other above 
mentioned relativistic correction has been reliably taken 
care of in the previous research on the atomic PNC prob- 
lem, see the references discussed above. The present 
paper demonstrates that all radiative corrections have 
recently been taken into account and calculated reli- 
ably. Thus, combining the radiative corrections with 
the previous theoretical data one can be absolutely cer- 
tain that every possible relativistic correction in many- 
electron atoms is accounted for. This done, the theory 
comes up with an important conclusion |48L l55| : the most 
accurate available experimental data for parity noncon- 
servation of Wood et al [JjJ fully support the standard 
model. 



II. VACUUM POLARIZATION AND RELATED 
PHENOMENA 

A. Variation of the electron wave function due to 
short-range potentials 

The effective Hamiltonian -ffpNC induced by the Z- 
boson exchange that describes the PNC part of the weak 
interaction between the atomic electrons and the nucleus 
can be presented in the following form (the relativistic 
units h = c = 1, e 2 = a are used, if not stated otherwise) 

ffpNC = ^G F Q w p(r)75 • (2.1) 

Here Gf and Qw are the Fermi constant and the nuclear 
weak charge, and p(r) is the nuclear density. The PNC 
amplitude describes the matrix element that mixes two 
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states of opposite parity for a valence atomic electron. 
According to Ea. (|2.1l) the electron wave function in this 
matrix element must be taken on the nucleus. Therefore 
the largest PNC amplitude arises due to the mixing of 
the s 1/2 andpx/2 partial waves (Vv, s ,i/2|#pnc|Vv, p ,i/2), 
here n, n! refer to the electron energy states, the electron 
wave functions are the Dirac spinors 



/ n y(r)O i/Al (n) 
ig n ij(r)n jIfl (n) 



(2.2) 



fnijif) and g n ij{ T ) are the large and small radial com- 
ponents of the spinor, f^^n) and fi .^- (n) = — (er • n) 

f2ji m (n) are the spherical spinors, and I + I = 2j. 

In applications the electron wave function is known in 
some reasonably good approximation (say, the Hartree- 
Fock-Dirac initial approximation plus contributions of 
the many-electron correlations). Inevitably there exist 
a number of perturbations that are not included in this 
approximation and can affect the electron wave function. 
The QED vacuum polarization is one of such perturba- 
tions, some others are considered below. In many cases 
the perturbation is localized in the vicinity of the nu- 
cleus, at distances r that are much smaller than the Bohr 
radius of the K-shell r < 1/ (aZm) . This is definitely 
the case for the QED vacuum polarization localized at 
r < 1/m. For this range of distances the behavior of the 
electron wave function is quite simple because it is gov- 
erned mostly by the pure Coulomb field of the nucleus. 
One can anticipate therefore that the influence of the per- 
turbation on the electron wave function can be presented 
in simple terms, directly via the perturbative potential. 
This anticipation proves true, there exists a simple way 
to account for the perturbation that was suggested in 
Ref.|3. 

Note that it suffices to find the relative variation of 
the electron wave function induced by a perturbation. 
Moreover, for small distances r <C l/(aZm) the relative 
variation does not depend significantly on r. Therefore 
it is sufficient to find the relative variations of the wave 
function at the origin, i. e. to find Sfnlj(0)/fnlj(ty an d 
5g n lj(0)/9nlj(0), where Sf n ij(r),Sg n ij(r) are the varia- 
tions of the large and small components of the wave func- 
tion due to the considered perturbation, and the relative 
variations at the origin are taken in the limit r — > 0. 
The limit is necessary because either the large or small 
components of the wave function vanish at the origin, see 
Ea. (|3.16|l below. Ref. JijjJ found for the relative varia- 
tions of the wave functions the following result 



Sf(0) Sg(0) m 



/(0) 5(0) 



^ / V(r)(a + kr)dr . (2.3) 



Here and below the indices nlj for the wave functions are 
suppressed to simplify notation, a is a parameter with 
length dimension and k is a dimensionless coefficient 



7 mc 
2k(2k - I) 

7(47 2 - 1) ' 



(2.4) 
(2.5) 



V(r) is the perturbative potential. Parameters 7, n are 
defined conventionally 7 = (1 — (Za) 2 ) 1 ^ 2 and k = 
K i + 1 )-i(i + 1 )- 1 / 4 =±(i+ 1 / 2 )- Relations (IP) . i2~4l 
arc presented in absolute units, to make them more con- 
venient for different applications below. The simple, clear 
formula 1)2. 3fl solves the problem formulated above, pre- 
senting a variation of the wave function at the origin 
in transparent terms, over the two first moments of the 
potential. It is instructive to consider Ea. H2.3fl in the 
nonrelativistic limit that reads 



h 2 21 + 1 



V(r) r dr 



(2.6) 



Here ip(r) is the Schroedinger nonrelativistic wave func- 
tion. Deriving Eq. (|2.6|) from (|2.3I) one uses the fact that 
according to (|2.4fl the parameter a goes to zero in the 
limit Za -> 0, while from 0110 nnds k ~> 2 /( 21 + 1 )- 

There is a simple short-cut derivation that leads to (|2.6|i . 
One expresses the variation of the Schroedinger electron 
wave function via the Green function G 



8ip(0) = I G(0,r')V(r')ij(r')d 3 r' , (2.7) 



approximating the Green function by the Green function 
for free motion 



G(r,r') 



Gf (r,r') 

2m 1 
~W 21 + 1 r l + 



(2.8) 



- Y lm (n)Y? m (n') , 



=-/ 



where Gf°^(r, r') is the Green function for the free mo- 
tion in the l-th partial wave, n = r/r, n' = r'/r', and 
the last identity in Ea. (|2.8(l takes into account the fact 
that the binding energy of the valence electron is negli- 
gible for short distances. Remembering that for nonrel- 
ativistic motion the wave function behaves like ipi( r ) (x 
Yi m (n) r , r — >• one immediately derives Eq. (|2.6I) from 
Ea. (1371). 

The point to be noted in this simple derivation is the 
dominance of the kinetic term in the nonrelativistic limit. 
For small distances the nonrelativistic kinetic energy be- 
haves as oc 1/r 2 , being larger than the potential energy. 
That is why Ea. p.Sfl holds in the nonrelativistic limit. In 
the relativistic approach the situation is different. The 
kinetic term behaves as oc 1/r, while the potential be- 
haves as Za/r. One has to expect henceforth that the 



5 



contribution of the Coulomb potential energy to the fi- 
nal answer is to be suppressed compared with the con- 
tribution of the kinetic term only by a factor Za. This 
is exactly what happens in Eq. H2.3f) . The second term 
oc kr on the right-hand side can be considered as a mere 
modification of the nonrelativistic result, while the first 
one oc a is the expected new term proportional to Za. 
Thus the structure of Eq. (|2.3f) is predetermined. The 
calculations of Ref. establish the coefficients in this 
relation as well as the dependence on higher powers of 
the Zen-expansion that appear through the parameter 7. 

The nonrelativistic Eq. if2.6f) shows that the parame- 
ter that governs the variation of the wave function is 
/ mV(r)rdr that is almost an obvious result valid for a 
variety of quantum mechanical problems |6l| . It is inter- 
esting that relativistic Ea. (f2.3f) shows that there exists 
another relativistic parameter J maV(r)dr. It is sup- 
pressed compared with the nonrelativistic one only by 
a factor Za that is not small for heavy atoms. This 
suppression can be well compensated for if the poten- 
tial increases at small distances which makes J maV (r)dr 
larger than m J V{r)rdr. In this case the new relativis- 
tic parameter becomes dominant, as the examples below 
demonstrate. The singular rise of the potential at small 
distances is typically smeared out by the finite nuclear 
size. Dealing with such singular potentials it is important 
to keep in mind that the first term on the right-hand side 
of Ea. H2.3fl is proportional to Z. Therefore it originates 
from those distances where the full strength of the electric 
field created by the nuclear charge is present, i. e. from 
the region outside the nuclear core r > ru- Inside the 
nucleus the Coulomb field created by the nuclear charge 
diminishes, making this region less important. In order to 
estimate its contribution, one can introduce for distances 
r < an effective nuclear charge Z — > Z e g = (r/r-^) 3 Z 
that diminishes with the nuclear radius taking into ac- 
count the decline of the electric field inside the nucleus. 
Correspondingly modifying the parameter a in Ea. (|2.3(l 



Here 6(ij} Si i^I-HpncIVv, 1/2) is * ne variation of the PNC 
matrix element due to variations of the wave functions 
created by the potential V(r). Deriving this result we 
take into account that essential parameters for the Si/ 2 
and pi/2 states are k s = —1, k p = 1, 7 S = 7 p = (1 — 
(Za) 2 ) 1 / 2 = 7 and define k in (|^TU|) to be k = k s = 

6/( 7 (4 7 2 -l)). 

It was assumed in Ea. (|2.3|l that the perturbative po- 
tential is a scalar. In applications it is necessary also to 
deal with more sophisticated, tensor-type perturbations. 
The electric and magnetic fields created by the nuclear 
dipole, quadrupole and higher moments give examples of 
such perturbations. It is interesting that Ea. (|2.3() can be 
generalized to describe the tensor-type perturbations in 
simple clear terms. Suppose that we have the valence 
electron described by the Dirac spinor ijjij (r) (the energy 
of the valence electron is so low that we can neglect it, 
suppressing the irrelevant index n). Suppose that there 
is some tensor perturbation V(r). When we consider the 
PNC amplitude we need to evaluate the wave function 
at the nucleus. There are only two partial waves that 
remain large at the origin, Sx/2 and P\/ 2 . This means 
that in the multipole expansion of the V (r)ijjij (r) term 
in the Dirac equation only these two partial waves con- 
tribute. The angular structure of these partial waves is 
predetermined, compare Ea. l|2.2|l . therefore we may write 



[F(r)^(r)k 1/2 



_ ( FL,l/2(r)tol/2,L,n( n ) 
iG i,l/2W fi l/2,i,/ n ) 



(2.11) 



where [V(r)ipij(r)] L1 / 2 means the admixture of the L 1 / 2 
state in the term V(r)ipij(r). Ea. (|2.11|l shows that the 
contribution of the term V(r)ipij(r) to the Dirac equation 
is expressed in terms of the scalar functions F L i/ 2 (r) 
and Gl i/ 2 (t) with L = 0,1. We can introduce now 
some effective potentials V£ ff (r) , Wf 3 (r) that satisfy the 
following conditions 



a c ff = — 



(2.9) 



one can apply this equation for the region r < m- We will 
see that this approximate procedure gives sensible results 
for the vacuum polarization, see Eas. H2.21fl - H2.23fl . 

Let us apply Eq. lj2.3fl to the PNC problem, in 
which one needs to calculate the PNC matrix element 
(^n, s ,i/2l#PNc|V>n',p,i/2)- Botn s i/2 and p 1/2 wave func- 
tions are influenced by the perturbation V(r). Therefore 
the relative variation of the PNC amplitude can be found 
simply by adding variations of S\j 2 and P\/ 2 states given 
in Eq. [j!Hjl 



JPNC, V 



^ S .l/2l#PNC#p, 1/2 
(lpa,l/2\HpNc\tpp,l/2. 

-wL v{r) 



(2.10) 



— kr 
3 



dr 



F LA/2 {r) = Vf(r)f Ltl/2 (r) , 
Gx,i/ 3 (r) = Wf(r)g L . 1/2 (r) 



(2.12) 
(2.13) 



Here /l,i/2(?"); 5i,i/2( r ) are the known large and small 
components of the nonperturbed wave function for the 
valence electron in the Li/ 2 partial wave. The proce- 
dure described shows that the term V(r)ipij(r) in the 
Dirac equation at small distances can be described by 
the effective potentials V£ s (r) and W^ s (r) (which can 
be different for the large and the small components for 
an arbitrary perturbation). Correspondingly, the varia- 
tion of the wave function 5ipij(r) at small distances can 
be expressed via its multipole expansion in which only 
the S1/2 and P1/2 states are retained 



#ii(r) 



L=0,1 



<5 /;j;i4/2( r )^l/2,L, A1 (n) 

a J!j;i,i/2( r ) !! i/2,i>) 



(2.14) 
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Here the functions <5/y ; l,i/2j $9ij;L,i/2 describe the ad- 
mixture of the L1/2 states to the function ipij(r) that 
arise due to the tensor-type perturbation. They are to 
be found from the Dirac equation. We see that the com- 
plicated multipole problem is reduced to the problem 
with the scalar potential, the only distinction is that the 
role of the scalar potential plays the effective potentials 
U£ ff (r), W£ s (r) defined above. However, this distinction 
makes no difference, we introduce the effective scalar po- 
tentials in the Dirac equation and solve it regardless of 
the nature of the potentials. We can rely on the discus- 
sion in to derive formulas similar to Ea. (|2.3|) . The 
result reads 



5fly,LA/2(0) 

/£,i/a(0) 

$9lj; £,1/2(0) 

#£,1/2(0) 



= —tfj V£ n (r)(a + kr)dr, (2.15) 

= — / Wf(r){a + kr)dr. (2.16) 
" Jo 



Remember that here Sfij-L, 1/2(0)1 $9lj; £,1/2(0) describe 
the admixture of £1/2, L = 0,1 states that arise in the 
wave function ipij (r) of the valence electron in the state 
with the quantum numbers Ij under the influence of the 
tensor-type perturbation V(r), while /z,,i/2(?"), 9l,i/2{ t ) 
describe the wave function of the valence electron in the 
state with the quantum numbers Li/2- The latter ones 
are the same as the ones used in Eas. (|2.12|) . <|2.13ll to de- 
fine the effective potentials. This makes our results pre- 
sented in Eas. I|2. 15|) . I|2. 16fl independent of the normaliza- 
tion of the wave functions /L,i/2( r )j 9l. i/2( r )- I n partic- 
ular, they are independent of the energy of the valence 
electron (as they should be since this energy is too low to 
produce any significant effects for small distances). The 
quantum numbers in Eas. (|2.4| ). (|2.5[) that define the pa- 
rameters a, k in the right-hand sides of Eas. (|2.15|l . (|2.16() 
correspond to the state £1/2, L — 0,1. We see that 
for an arbitrary tensor-type potential the result remains 
very simple. The only modification, compared with the 
simplest scalar case, is that one needs first to define the 
effective potentials using Eqs. (|2.12jl . l|2.13[) . 

Summarizing, we described in this section simple use- 
ful procedures that express the variation of the valence 
electron wave function at the origin directly in terms of 
the perturbative potential. 



B. QED vacuum polarization 



The vacuum polarization in the lowest order of pertur- 
bation theory is described by the Uehling potential Vyp 

M 



Vy P (r) 



no 



2a ( Ze 2 
3ir \ r 



c nc)dc, 



±\v?E± 
x 2 ) e ' 



(2.17) 



(2.18) 



Importantly, the Uehling potential l|2.17|l is singular at 
the origin, its asymptotic for mr <C 1 reads 



VVp(r) 



2a 
3^ 



Ze 2 
r 



In- 



C 



(2.19) 



the singularity is smeared out by the nuclear finite size 
r N . In Ea. (f2~T%f) C = 0.577 ... is the Euler constant. The 
function lnrnr arises due to conventional scaling of the 
QED coupling constant e 2 that manifests itself at short 
distances. Alternatively, this fact is referred to as the 
unscreening of the nuclear charge at small distances. For 
large distances mr 3> 1 the Uehling potential exponen- 
tially decreases 



V V p(r) = - 



Ze 2 \ a exp(— 2mr) 



r J 47T 1 / 2 (mr) 3 / 2 



(2.20) 



In order to find the influence of the vacuum polariza- 
tion on the PNC amplitude we substitute Ea. (|2.17|l into 
Ea. (|2.10|l . Integration over the variable r is performed 
analytically resulting in the following expression for the 
relative correction to the PNC amplitude 



where 



5 (vp, 



5 (VP, 



5 (VP) _ £(VP,1) + J(VP,2) 



Za 2 



(2.21) 



4 4 7 2 - 1 



(2.22) 



Za 
"27 



( VVp^n) n* ) 



(2.23) 



Here the function E\ (x) is defined conventionally 



£i(z) 



f°° dt 
J exp(-xt) — . (2.24) 



Ea. (|2.22|l . which was first found in gives the main 
contribution |77|. Deriving it one restricts the integra- 
tion of the singular term J maVyp in Eq. H2.10fl to the 
region outside the nucleus r > r^. The large Coulomb 
field of the nucleus makes this region important. Inside 
the nucleus r < tn the electric Coulomb field diminishes, 
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suppressing the contribution of this region. In order to es- 
timate the contribution of small distances r < r n we use 
a modification of Eq. 1)2.10(1 introducing into Eq. Il2.3|) the 
effective parameter a defined in Ea. l|2.9|l . The resulting 
integral in Eq. (|2.3I) is saturated in the nuclear interior 
region in the vicinity of the nuclear surface, where the 
electric field remains large. This fact allows one to take 
the values of 7 and Vyp(r) directly at the nuclear surface. 
The described procedure results in estimate Ea. H2.23|) . 

Eas. (|2.21fl - H2.23[l give the weak interaction matrix el- 
ement for an arbitrary atom in a transparent analytical 
form without fitting parameters. Numerical results are 
easily obtained by a straightforward one-dimensional in- 
tegration in Ea. (|2.22J) . One only needs to specify the 
nuclear size that can be taken conventionally as tm — 
1.1 A 1 / 3 fm, where A is the atomic number |78|, Alterna- 
tively, the right-hand side of Eq. (|2.21|) can be calculated 
using an expansion in powers of mr^ <C 1 that reads |43| 



5 (VP) = « 
7 

h {Zaf 



4(4 7 2 
1 



1) 



Za- 



(2.25) 



In 



mr/v 



0.759 



where the omitted terms are 0{mr^) |79j. The error 
of simplification of Ea. (|2.10|l to (|2.25(l remains < 7% 
for heaviest atoms, decreasing for lighter atoms. A 
very interesting feature revealed by Ea. H2.25fl is the 
second power of the logarithmic factor in the term ~ 
Z 2 a 3 In 2 (mr^). This strong, double-logarithm enhance- 
ment was first found by Milstein and Sushkov 0t The 
derivation presented above makes the origin of this ef- 
fect clear. The first power of the logarithmic func- 
tion In mr originates from the unscreening of the nu- 
clear charge at small separations (|2.19|l , while the second 
one appears from the relativistic parameter J V{r)dr oc 
J drhi(mr)/r oc In 2 1/mr^ in Ea. (|2.10|l . 

In Fig. n we present the relative corrections due 
to the QED vacuum polarization <5( VP ) to the PNC 
matrix element calculated with the help of Eas. i|2.21|) - 
(|2.23|) . The corrections are positive and rapidly grow 
for heavy atoms. For the most interesting case of the 
133 Cs atom Ea. (|^2T|) give s <5< VP ) = 0.40%. For 205 T1, 
208 Pb, and 209 Bi Ea. l(^T|) predicts corrections <5 (VP ) = 
0.93, 0.96, 0.99% respectively. For all four atoms the 
nuclear interior region gives 10 % of the total correction, 
£(vp : 2) _ o.lt)( vp ), which means that to make results 
accurate numerically this region needs to be taken care 
of. 

Compare these results with other results reported re- 
cently. Johnson, Bednyakov and Soff 0] calculated for 
the first time the correction due to the vacuum polar- 
ization for the 6s — 7s PNC amplitude in 133 Cs and 
found it to be 0.4%. The result of [4l| includes, along 
with the variation of the weak matrix element, variations 
of the dipole matrix element and the corresponding en- 
ergy denominators that, combined together, describe the 



S1/2 — Si /a amplitude measured experimentally. Dzuba 
et al [36j confirmed this result, calculating the correc- 
tion 0.41%, and supplied more details providing sepa- 
rate variations for all three quantities mentioned above. 
They found that variations of the dipole matrix elements 
and the energy denominators induced by the Lamb shift 
corrections, being not small, compensate each other al- 
most completely. Thus the variation of the PNC am- 
plitude proves to be equal to the variation of the PNC 
matrix element. The numerical result for thallium found 
in 66] gives a correction 0.94%. to be compared with 
0.93% mentioned above. Calculations of Milstein et al in 
Ref. 0] were restricted by the logarithmic accuracy that 
was further improved by using some constant as a fitting 
parameter to reproduce 0.4% for 133 Cs, in line with [4l| . 
For Tl Ref.jH finds 0.9%. This is slight ly less than the 
above mentioned prediction of Ea. H2.21|) (0.93%). This 
(insignificant) discrepancy probably arises because the 
procedure adopted in [I^] does not account for the in- 
crease of the contribution of the nuclear interior region 
with the nuclear charge. We conclude that all reported 
calculations for the vacuum polarization correction are 
in very good agreement, distinctions between results de- 
rived by different methods and by different groups are all 
below 5%. 

Let us discuss now possible corrections to the results 
discussed above. One of them originates from a modifica- 
tion of the vacuum polarization due to the finite nuclear 
size. The corresponding generalization of the Uehling 
potential was given by Fullerton and Rinker |63j. In or- 
der to estimate its contribution we need to remember 
that the main oc In 2 mr^ term in Eq. (|2.25|l originates 
from the integration well outside of the nuclear radius 
where the Uehling potential is very close to the one of 
Fullerton and Rinker. Therefore this main part of the 
correction is not changed. Modifications could arise only 
for the integration in the vicinity of the nuclear surface, 
but this region by itself gives only a fraction (~ 10%) of 
the total correction, while the relative difference of the 
Uehling potential and the Fullerton and Rinker potential 
on the nuclear surface is small. For all atoms it remains 
< 5%, as we found by direct numerical comparison of the 
Uehling potential with the Fullerton and Rinker poten- 
tial. We conclude that the modification of the vacuum 
polarization due to the finite nuclear size is irrelevant. 

Another possible modification can originate from the 
QED higher-order corrections that contribute to the vac- 
uum polarization. One of them stems from the fact 
that the Landau pole requires that the linear in In mr 
term in Eq. (|2.19|l be followed by higher-order logarithmic 
terms. However, the corresponding series runs in powers 
of Za{a\nmr) n , n — 1, 2 ... , giving small contributions 
for higher order terms n > 1. 

Consider now a modification of the Uehling poten- 
tial that is due to the fact that the Coulomb field of 
the nucleus disturbs the virtual electron-positron pair. 
This effect is conveniently described with the help of 
the Wichmann-Kroll potential [64| that, compared to 
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Ea. (|2.17|l . describes the next term in the Za-expansion. 
For estimates it is sufficient to take the Wichmann- 
Kroll potential at small distances, where its asymptotic 
mr <C 1 reads 65): 



2a 

V WK (r) - 0.092— [Zaf 
on 



Ze 2 
r 



(2.26) 



Using Eq. (|2.10|l one immediately finds from Eq. I|2.26|l the 

ratio of the correction induced by the Wichmann-Kroll 
potential to the Uehling correction 



x(VP) 
°WK 

§ (VP) 

°u 



= -0.184- 



(Zaf 



In l/(mr N ' 



(2.27) 



as was found in [43. For 133 Cs this ratio is about -0.007, 
demonstrating that the correction due to the Wichmann- 
Kroll potential is very small. Thus all higher order QED 
terms in the vacuum polarization are negligible. 

Alongside the vacuum polarization there exists also an- 
other polarization effect, namely the polarization of the 
nucleus. One can estimate its influence on the electron 
wave functions and the PNC matrix element with the 
help of the conventional polarization potential V^p that 
is created by the nucleus 



e 2 a d 



2r 4 



(2.28) 



is of the order of ~ 0.1Z 2//3 a 2 jgjj, which numerically is 
close to Eq. lESP . 

Summarizing, we demonstrated above that the QED 
vacuum polarization gives a large positive correction 
(~ 1%) to the PNC amplitude. Contributions of higher 
order QED corrections are much smaller. This means 
that the influence of the nuclear Coulomb potential on 
the virtual electron-positron pair is not pronounced, pos- 
sibly due to opposite signs of the virtual particles. We 
will see below that for a single electron the situation is 
different, the Coulomb field does play an important role 
for the self-energy correction, see Eq. l|3.29l) . The influ- 
ence of the finite nuclear size on the vacuum polarization 
was found to be nonessential, the polarization can be well 
described by the simplest Uehling potential. Also negli- 
gible is the polarization of the nucleus. We presented 
above in some detail analytical methods of calculation, 
paying less attention to purely numerical methods that 
have played a very important role in the problem, see 
[36I IlO . l6t| . Omissions of details in our presentation of 
numerical methods is justified by the fact that the vac- 
uum polarization is described by a sufficiently simple po- 
tential, which can be incorporated in the initial Hartree- 
Fock-Dirac approximation for the atomic electrons, and 
then dealt with by conventional methods of many-body 
theory. Thus the numerical approach follows the main- 
stream of atomic structure calculations discussed in some 
detail elsewhere, see for example [32>|33|]. Overall, the po- 
larization is a well understood, robust effect that is well 
described by Eas.l|2"2T )) -l|!£55 )) . 



Here ad is the static nuclear dipole polarizability. It 
can be estimated roughly assuming that the total os- 
cillator strength for the nuclear dipole excitations is lo- 
cated at the energy of nuclear giant dipole resonance fid- 
Then, using the conventional sum rule for the oscillator 
strengths one finds 



a d 



Ze 2 

Tip 



(2.29) 



where m p is the proton mass. Since the potential 
increases at small separations the main contribution in 
Ea. (|2.10|l is given by the first term in Ea. (|2.10|) which 
should be integrated outside the nuclear radius r > r N , 
as was explained above. Having this in mind, one substi- 
tutes (|2.28(l in (|2.10() finding an estimate for the influence 
of the nuclear polarization on the PNC matrix element 



III. SELF-ENERGY AND VERTEX 
CORRECTIONS 

A. Chiral Invariance and QED corrections 

Let us show, following Ref.^^|, that there exists an 
interesting and useful relation that gives the QED ra- 
diative corrections to the PNC matrix element via simi- 
lar radiative corrections to another quantity, the energy 
shifts of the atomic electron induced by the finite nuclear 
size (FNS). This relation can be written as 



- (^FNS, s 



), 



(3.1) 



where <5pnc. sp is the relative radiative correction to the 
PNC matrix element between the Sx/2 and P1/2 orbitals 



<j(NP) _ ± 



3 7 to p Qjj r% 



(2.30) 



For the 133 Cs atom we take f2d — 15 MeV obtaining 
£(np) ^ o.01% that shows that the nuclear polarization 
is negligible. This conclusion agrees with estimates of 
Milstein and Sushkov 01 > wno proposed that the effect 



3pNC,sp 



(V> s ,l/2|gpNc|^,l/2) rad 
(V' s ,l/2|-ffpNc|V , p, l/2> 



(3.2) 



The energy difference between the considered opposite 
parity orbitals ~ 1 eV is much lower than a typical exci- 
tation energy ~ m = 0.5 MeV that governs the radiative 
corrections. We can therefore neglect this difference as- 
suming that E s> i/2 — E p i/2- This assumption makes the 



9 



correction (ip St \/2\Hpnc\^ p , i/2) rad gauge invariant. To 
make this argument more transparent, one can use the 
Coulomb approximation for the atomic field and consider 
the degenerate Coulomb nsi/ 2 and np 1 / 2 levels. 

The problem of gauge invariance can be viewed from 
another perspective. One can consider the full PNC am- 
plitude. Alongside the PNC matrix element it includes 
also the conventional El amplitude and the correspond- 
ing energy denominator. The PNC amplitude describes 
the events on the mass shell, the photon frequency in 
the El transition equals the excitation energy of the 
atom. Therefore the amplitude is definitely gauge invari- 
ant. Considering the radiative corrections in the lowest 
order, one can separate them into corrections to the PNC 
matrix element, the Lamb-shift corrections to the ener- 
gies, and the corrections to the El transition. Their sum 
is definitely gauge invariant. The point is that the lat- 
ter two corrections are known to compensate each other 
almost completely |36| . The only remaining corrections 
are the ones which are associated with the PNC matrix 
element. Thus these latter corrections are gauge invari- 
ant. To make this discussion complete, let us mention 
also that in the perturbation theory approach outlined 
in the next section llll Bl the gauge invariance of the PNC 
matrix element is self-evident. There is no doubt, there- 
fore, that one can assume that the radiative corrections 
to the PNC matrix element are gauge invariant. 

The quantities <5fns,s and 5fns,p m Ea. (|3.2l) are the 
relative radiative corrections to the FNS energy shifts 
-S'FNS.s; -^fns,pi for the chosen Si/ 2 and p\j 2 electron 
states 

^rad 

^FNS, I — ? / = 0,1. (3.3) 

-C'FNS.i 

The FNS energy shifts can be presented as matrix ele- 
ments of the potential VfnsM that describes the devia- 
tion of the nuclear potential from the pure Coulomb one 
inside the nucleus due to the spread of the nuclear charge 

£fNS,( = 1/2|T / FNS|V'2, l/2> , i = 0, 1. (3.4) 

Equality i|3.1[l may be established for the sum of all QED 
radiative corrections (~ Za 2 f(Za)) 1 or specified for any 
gauge invariant class of them. We will discuss first the 
self-energy and vertex corrections, restricting our consid- 
eration to the lowest order of perturbation theory pre- 
sented by the Feynman diagrams in Fig[3 calling them 
the e-line corrections. Later we will switch our atten- 
tion to the vacuum polarization and demonstrate that it 
satisfies (I3.1|) as well. 

The intermediate electron states in diagrams of Fig. [21 
are described by the propagator G = (j^p^+joU—m)^ 1 , 
with p^ — (e, — iV), where e — m + E m is the rel- 
ativistic electron energy, U — U(r) is the atomic poten- 
tial that includes the potential created by the nucleus 
with FNS and, generally speaking, the screening poten- 
tial of atomic electrons. The latter one is not important 



at small distances, but it still can be accounted for by 
the formalism. For the relativistic propagator we use the 
same symbol G that was applied previously for the non- 
relativistic Green function in Ea. (|2.7|) . which should not 
produce confusion. At small distances the exchange po- 
tential is negligible, therefore the potential is local. The 
external legs of the diagrams describe the wave functions 
"0s,i/2( r ) and ^ Pj i/2(r) for the considered Si/ 2 and pi/ 2 
levels. Let us note first that in the region of short dis- 
tances the following relations hold 

75 G = -G 75 , (3.5) 
VV/2O") = c75V , P ,i/2( r ) ■ (3-6) 

They follow from the chiral invariance that manifests it- 
self at small distances. For r <C 1/m the mass term in 
the Dirac equation becomes smaller than the kinetic term 
(oc l/r). For even smaller distances r < Za/m the po- 
tential (Za/r) also exceeds the mass term. Therefore for 
small distances we can neglect the mass m 0, using 
the resulting chiral invariance for our advantages. Since 
Eas. l|3.5(l . l|3.6(l will be important for us, we examine them 
in some detail. In order to prove Ea. l|3.5|l one neglects 
the mass term in the electron propagator assuming that 

G- 1 = 7^ +loU . (3.7) 

Then, the identity 757^ = —7^,75 leads to G~ 1 7 A , = 
— 7 M G _1 that results in l|3.5|l . In order to verify Eq. i|3.6|) 
let us write down the conventional Dirac equations for 
the large / and small g components of the electron wave 
function [|2.2[l 

/' + — / -(e + m-U)g = , (3.8) 
r 

9' + —.9 + (e-m-U)f = . (3.9) 
r 

Here we omitted the indexes njl, distinguishing the states 
by the parameter k that takes the values —1 and 1 for the 
S1/2 and p\/ 2 states respectively. For small separations 
one can neglect in Eas. (|3.8| . (|3.9|l the mass term m. Af- 
ter that we observe that the equations become invariant 
under a substitution 

/ - 9, (3-10) 
9 - -/• 

Remember now that according to Ea. (|2.2|) the s 1 / 2 and 
Pi/ 2 wave functions have the form 

p,/ V l 9p,l/2U J 
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Here the two component spinor u = u a , a = 1, 2 depends 
on the projection of the total momentum /i, u a = 5 a i 
for fi = 1/2, and u a = 5 a ,2 for [i — —1/2. It follows from 
Eci. (|3~TT|> that 

(-i/„,f;^".„)„) . <" 2 » 

where we took into account that in the standard repre- 
sentation considered for the Dirac matrixes 

75=(;j)- (3-13) 

From Ea. [j3.12|) . [j3.11j) one observes that Ea. [j3.6j> can be 

presented as 

/s,l/2 = c 5p,l/2 > (3-14) 

08,1/2 = -c/p,l/2 ■ (3.15) 

We see that these latter relations coincide with Ea. (|3.10|l . 
The first substitution k — ► — K in (|3.10() exchanges equa- 
tions for the S1/2 a nd px/ 2 states, while the two others 
are identical to Eas. lj3.14j) and l|3.15|l respectively. Thus 
the verified above Eas. (|3.10|l prove validity of Ea. lj3.6jl . 

At this point one may recall a subtlety, a symmetry 
of the Dirac equation revealed by Ecis. H3.10|) is, generally 
speaking, not sufficient. One needs to verify also that 
the boundary conditions satisfy same symmetry as well. 
However, one can argue that the boundary conditions 
arise as a physical complement to the Dirac equation and, 
therefore, the symmetry of the equation should automat- 
ically result in the symmetry of the boundary conditions. 
To verify this argument one can recall that for the finite 
nuclear charge distribution the boundary conditions at 
the origin read simply 

f /.,i/a(0) = const , f / p , 1/2 (0) = , 
1 03,1/2(0) = , \ fltp, 1/2(0) = const , 

(3.16) 

where const 7^ 0. Eas. (|3.16|l explicitly demonstrate that 
indeed, the boundary conditions are invariant under the 
substitution Ea. (|3.10|l . This completes our verification of 
Eci. (j3.6|) . The constant c in 1)3. 6|) . which depends on the 
normalization conditions, will be irrelevant for us because 
we are interested in relative quantities that appear in 
Ea. (|3.1|l where this constant will be canceled out. 

It is instructive to look at numerical verifications of 
Ea. H3.6|l shown in Fig|3] Supporting the proof given 
above, they provide also details that are cumbersome 
for analytic treatment. FigOl (a) shows that Eo. (|3.14jl . 
which according to the discussion above holds at small 
distances r -C 1/m, remains approximately valid at suf- 
ficiently large distances, up to (and even beyond) the 
Compton radius r < 1/ra. In contrast, Eo. (|3.15jl is ap- 
plicable only at small r <C 1/m as shown in Fig. [3] (b), 



being violated at r ~ 1/m, as one can see from Fig. Ola). 
We can remember, however, that the dominant contribu- 
tion to the PNC matrix element arises from the / Sj i/2 and 
<?p,i/2 functions, which are both large for r < 1/m, much 
larger than the other two functions g S: i/2, fp, 1/2- Since 
these large functions satisfy Eq. (|3.14|) even at sufficiently 
large distances r < 1/m, one can also rely in this region 
on Ea. H3.6|) . We will exploit this opportunity at some 
later stage below, but our prime interest is related to the 
distances as small as the nuclear radius r < r-^ <C Za/m, 
where we verified beyond doubt that Eos. (|3.5j l. (|3.6jl hold, 
see FigGU(b). 

Consider the diagrams (b),(c) in Fig. [21 that describe 
the self-energy correction for the PNC amplitude. The 
thick lines in this figure represent the electron propaga- 
tion in the atomic field, the dot vertex is the PNC inter- 
action described by Eo. (j2.1jl . and the dashed line shows 
the photon. Let us specify that the left and right legs 
of these diagrams describe the S1/2 and P1/2 states re- 
spectively. The PNC interaction Ijz.ljl is located at small 
distances tn <C 1/m. Let us use this fact and, relying 
on Eos. (|3.5j ) (|3.6|) . transform the PNC matrix element. 
To simplify the notation it is convenient to present the 
Hamiltonian of the weak interaction Eq. (j2.1j) in the form 

-HpNC = V PN c(r)j5 , (3-17) 

where VpncM = Gf Qw p(r) / (2^/2) is a scalar factor 
in the PNC interaction. Consider first the diagram (b). 
The f/'p,i/2( r ) function in this diagram enters the PNC 
interaction that makes its argument r localized inside the 
nucleus, r < rj$. This allows us to use Eo. H3.6jl deriving 

(^,i/2l75%cGS|^ 1/2 ) (3.18) 
= {tps, 1/2I Vpnc G E|V> S , 1/2) ■ 

Here and below we use the symbol of identity (=) that 
simply means dropping the irrelevant constant c intro- 
duced in Eo. (j3.5|) (that is canceled out when the rela- 
tive contribution of the correction is considered anyway). 
Ea. (j3.18jl uses conventional notation, S and G are the 
self-energy operator and the electron propagator men- 
tioned above. The right-hand side of this relation shows 
that we can look at the diagram (b) as the one that de- 
scribes the diagonal Si/ 2 — Si/ 2 transition induced by an 
effective scalar potential Vpnc- Consider now the same 
diagram (b), but for another physical quantity, for the 
FNS energy shift for the Si/ 2 level. Its contribution to 
the energy reads 

(^.,i/ 2 |VFNsG£|^ fl/9 > • (3.19) 

Eos. lj3.18l) . lj3.19jl give the contribution of one and the 

same diagram (b) to different quantities, namely the 
PNC amplitude and FNS energy shifts. Comparing them 
one observes their similarity. In both cases the diagram 
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describes the s\/2 — S1/2 transitions that are induced by 
the short-range scalar potentials. The only distinction 
comes from the fact that the potentials VpNc( r ) an d 
VfnsM possess different profiles. 

This latter difference, however, proves to be irrelevant. 
In order to see this we need to discuss the integral with a 
short-range potential that arises in the diagram. Its inte- 
grand includes the potential (i. e. VpncM f° r the PNC 
amplitude and VfnsM f° r the FNS energy shift), the 
electron wave function ips,i/2{ r ) & n d the electron prop- 
agator G(r,r'). The wave function and the propagator 
inside the nucleus are almost constants, see discussion of 
this fact below. Taking them out of the integration we 
observe that the considered sophisticated integral is sim- 
plified, being reduced to the integral over the potential 
itself J VpNc{r)dr. It is easy to see that the same integral 
arises in the non-perturbed matrix element. Therefore in 
the relative contribution of the correction this integral is 
canceled out, and the result does not depend on the po- 
tential at all. This conclusion is valid for any short-range 
potential. Thus the relative contribution of the diagram 
(b) proves to be one and the same for the PNC and the 
FNS problems. This means that the relative contribution 
of the diagram (b) to the PNC matrix element equals half 
of its relative contribution to the energy shift of the Si m 
state due to FNS. The mentioned factor 1/2 takes into 
account the fact that the energy shift is induced by the 
two diagrams (b) and (c) in Fig. |2]that are identical for 
the FNS problem. 

Before proceeding further, let us discuss the behavior 
of the electron wave functions inside the nucleus which 
was exploited in the above derivation. Note that inside 
the nucleus the large functions f St \/t(r) and g Pt i/2(r) ex- 
hibit smooth, almost constant behavior, as demonstrates 
Fig. El The reason for this effect stems from the fact that 
inside the nucleus the electric field rapidly diminishes 
with the radius. The functions g s ,i/2( r ) an( i /p,i/2( r )i 
show some variation, but it remains very mild. Assessing 
consequences of this smooth variation one can recall that, 
firstly, the integration volume r dr enhances the contri- 
bution of the nuclear surface region where these functions 
are very smooth indeed, as shows Fig. 01(b). Secondly, 
their role is not prominent anyway, simply because they 
are smaller than the functions f s ,i/2( r ) 9p. i/2( r )- There- 
fore they can be well approximated by constants. Thus 
both s\/2 and P1/2 states are well described by constant- 
type functions inside the nucleus. This fact was used 
above for the Si/ 2 state, the argument is repeated be- 
low for the Pi/2 state. We also presumed above that the 
propagator G(r,r') does not depend on the variable r, 
when r < rpj- This statement follows from the simple 
fact that the coordinate r' gives the location of the point 
where the radiation of the photon takes place, therefore 
r' ~ 1/m 3> r. The smallness of r compared with r' 
justifies its neglect in the propagator. 

Examining the diagram Fig[2] (c) we use the same ar- 
guments as we exploited above for the (b) diagram prov- 
ing that its relative contribution to the PNC amplitude 



equals half of its relative contribution to the energy shift 
due to FNS of the pi/ 2 state. Combining our results for 
the diagrams (b) and (c) we conclude that the self-energy 
corrections (superscript SE) described by these diagrams 
comply with Eq. H3.lfl 

r(SE) _W 5 (SE) -(SB) \ ,„ 2Q s 

°PNC, sp — 2 \ °FNS,s ^ u FNS,p J • 

The given derivation does not use any specific gauge con- 
dition, therefore it is gauge invariant. 

Consider now the vertex correction (a) in Fig. [21 An 
analytical expression for it can be written in the conven- 
tional form 

-iaJ-^drdv'dr D^(q) (3.21) 
xe"'( r - r ')$ F (r,r',ro) , 

where 

$ MIy (r,r',r ) = ^ p ,i/ 2 (r)7 AI G(r, r ; e + w) 

x VpNc(?'o)7o75G(ro,r / ;e + ^) 

x 7^ s ,i/ 2 (r') • (3.22) 

Here D^ v (q) is the photon propagator, = (u;,q) is 
the four-vector of the photon momentum. The range of 
distances r ra d where the radiation processes take place is 
of the order of the Compton radius, r ra d ~ m . This 
ensures that r, r' ^S> r$ because the PNC interaction is 
localized at ro < r^. We can simplify therefore Ea. H3.21[) 
factorizing the interaction over ro 

J dr 3v(r, r',r ) = / VpNc(^o)dr 

xii/ 2 (r)7,G(r, 0;e + W ) (3.23) 
x 7o 75 G(0, r'; e + u) 7^ ^ s ,i/ 2 (r') , 

where we put ro = in the arguments of the propaga- 
tors. We observe that the potential arises only in a factor 
/ VpNc( r o) dro that is separated from all other elements 
of the diagram. The same factor appears in the main 
term of the PNC amplitude. Therefore the relative con- 
tribution of the vertex correction does not depend on the 
potential at all. Remember that the same property is 
exhibited by the self-energy correction discussed above. 

The variables r, r describe the radiation processes and 
therefore are expected to be of the order of the Compton 
radius r, r' ~ 1/m. Let us consider the implications that 
arise if one presumes that the region of small separations 

r,r'<— (3.24) 
m 

gives the dominant contribution to the integral in 
Ea. (|3.21|l . We can engage in this region Eas. (|3.5|l . (|3.t)|l . 

The first of them allows one to exchange the 75 matrix in 



12 



Eq. (|3.23|l with the propagator G. One can also transpose 
it with the vertexes 757^ = —7^,75- Using several such 
transpositions we can bring the 75 matrix from its central 
position in expression (|3.23[) to the side, where it hits the 
external wave function. After that we use Ea. (|3.6|) for 
this function. Repeating the procedure twice, one time 
shifting 75 to the left, and another one to the right we 
derive 



Vi p (r)7 M G(r, 0) 7o 7 5 G(0, r')7^(r') 

= fo(r) 7 p G(r, 0) 7o G(0, r') 7 „ ^ s (r') (3.25) 

= ^ p (r) 7m G(r, 0) 7o G(0, r') 7, ^ p (r') . (3.26) 

To simplify notation we dropped here the total mo- 
mentum 1/2 from the indexes in the wave functions 
(i. e. ip p = ip p ,i/2 etc), and omitted the energy vari- 
able in the propagators (this variable equals e + ui for 
all propagators). The point of these transformation is 
that the 75 matrix disappears from the right-hand sides 
of Eqs. (|3.25[l . (|3.2t)[) . Substituting these two expressions 
back in Eas. H3.23|) . (|3.21|l we find that the relative con- 
tribution of the vertex diagram (a) to the PNC matrix 
element can be expressed either as a matrix element of 
the S1/2 — S1/2 transition, or the matrix element of the 
P1/2 ~ P1/2 transition. We recognize in these two matrix 
elements the relative contributions to the energy shifts 
due to FNS for the S1/2 and P1/2 states respectively. This 
latter conclusion takes advantage of the fact discussed 
above that the profile of the short-range potential does 
not influence the relative correction, one can conveniently 
choose the potential to be Vpnc or Vfns • This discussion 
demonstrates that the vertex (superscript V) corrections 
to the PNC amplitude and to the energy shifts due to 
FNS prove to be approximately equal 



r(V) ^ AV) ^ ,(V) 
°PNC,sp — °FNS,s — °FNS, 



(3.27) 



Combining Ecis. H3.20|) . H3.27fl we nn< ^ that the total e-line 
corrections, i. e. the self-energy plus the vertex correc- 
tion, satisfy Eq. l|3.1|l . 

The derivation presented above exploits Eq. (|3.24|l pre- 
suming that the small distances are dominant in the ver- 
tex correction. There are reasons indicating that this as- 
sumption does not put restrictions on the derived result 
(|3.1() . As a first attempt one can try to exploit the be- 
havior of the wave functions / s l / 2 , 9 P .i/2 - Fig- EDshows 
that they approximately satisfy Eq. l|3.6[l even at suffi- 
ciently large separations, of the order of the Compton 
radius r < 1/m. This shows that Ea. H3.25[) remains ap- 
proximately applicable even if the integration region in- 
cludes sufficiently large separations r ~ 1/m. This, in 
turn, demonstrates that the first equality in Ea. H3.27fl . i. 
e. Sp^ c sp — <5p^g s also remains approximately valid. 
This is a positive result. However, this fact alone is not 
sufficient. The second equality in Ea. H3.27[) proves to be 
much more sensitive to the essential integration region 



because the wave functions f p ,i/2, 9s, 1/2 fail to satisfy 
Ea. (|3.6|) at r ~ 1/m. We need therefore to be more care- 
ful. 

Let us recall at this stage that the range of distances 
T-cad where the radiation processes take place depends 
on the chosen gauge for the electromagnetic field. It 
is well known that the usual length-gauge favors larger 
distances, while the velocity and, especially, acceleration 
forms, make the smaller distances contribute more to the 
radiation process. These known examples show that the 
region that contributes to the radiation process is not 
gauge invariant. This is important. If the radiation re- 
gion is not gauge invariant, then we are able to choose 
a gauge in which this region is located close to the nu- 
cleus, inside the zone defined by Ea. (|3.24|l . In this gauge 
the above derivation is valid and Ea. 13.27f) for the ver- 
tex is correct. Note that this argument does not rely on 
an explicit form for the necessary gauge. It suffices to 
acknowledge only that the radiation region is not gauge 
invariant. Summarizing, there exists a gauge (more ac- 
curately, a family of gauges) in which Ea. H3.27[) holds. 

Recall now that Eq. 1)3.20(1 for the self-energy remains 
valid in any gauge, in particular in the one discussed 
above. In this gauge we can combine together the self- 
energy and vertex corrections (|3.2(J|) . I|3.27|I proving the 
validity of their sum Eq. l|3.1|) . The three quantities in 
this latter relation are all well-defined, gauge invariant 
physical observables. The fact that Eq. 13.lfl is derived in 
one particular gauge means that it is valid in any gauge, 
as was proposed in |48|. 

Let us summarize the above discussion. We argued 
that the calculations can be organized in such a way 
that all important events happen at small distances. 
In this region the chiral invariance holds, resulting in 
Eqs. H3.5|) . l|3.6l) and, consequently, in Eq. l|3 1|) . Thus this 
latter equality expresses the fundamental and simple fact. 
At small, nuclear distances the chiral invariance governs 
the problem. Consequently, Eo. (|3.1|) that follows from 
it, can be called the chiral invariance identity. 

Let us estimate the accuracy of Ea. (|3.1ll . The above 
derivation used the gauge in which the radiation pro- 
cesses for the diagram (a) take place mostly at small dis- 
tances r ra d < m . They should also take place outside 
the nucleus tn < r m d, as is necessary to justify our pre- 
sumption that the shape of the potential inside the nu- 
cleus is irrelevant. From the last two inequalities we find 
that the derivation relies on a parameter £ = mr-^ ~ 0.01. 
This determines the magnitude of the error (few per cent) 

of Ea. cm. 

Similarly one considers the contribution of the QED 
vacuum polarization. Eqs. (|2.3| . (|2.4|l and (|2.5() present 
explicit variations for Si/ 2 and P1/2 wave functions at 
the origin induced by the vacuum polarization. Using 
these wave functions to calculate corrections to the PNC 
matrix element and FNS energy shifts, one immediately 
finds that the chiral invariance identity l|3.1|l holds for 
the vacuum polarization as well. This fact is in line with 
Eo. (|3.2U|l for the self-energy corrections. The proof of 
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(|3.1(l for the vacuum polarization stops here because one 
need not care about the complicated vertex corrections. 

Note that we do not consider above the radiative cor- 
rections of the order ~ a/7r which appear in the plane 
wave approximation. These contributions have been in- 
cluded into the radiative corrections to the weak charge 
Qw (and the renormalization of the charge and electron 
mass in the case of FNS energy shifts). Correspond- 
ingly, we subtract the contribution of the plane waves 
from Ea. (|3.1(l . considering only the part of the correc- 
tions that depends on the atomic potential ~ Za 2 f(Za). 
For heavy atoms this subtlety is insignificant numerically 
because the considered Z-dependent part of the correc- 
tion is bigger than the omitted Z-independent one, as we 
will see below. 

The chiral invariance identity l|3.1|l dives the e-line cor- 
rections to the PNC matrix element, which are difficult to 
calculate, in terms of the corrections to the FNS energy 
shifts that have been well-studied both numerically, by 
Johnson and Soff 0] , Blundell , Cheng et al I52J and 
Lindgren et al , and analytically, by Pachucki |67j ] and 
Eides and Grotch and Eides et al [j^ . Ref. [&2 presents 
the e-line radiative corrections to the FNS energy shifts 
for lsj/2, 2si/2 and 2p^/ 2 levels in hydrogenlike ions with 
atomic charges Z = 60, 70, 80, 90. Eq. (|3.1|l contains 
relative corrections, therefore we needed to calculate the 
FNS energy shifts -Efns ■ This was done in by solving 
the Dirac equation with the conventional Fermi-type nu- 
clear distribution p(r) = pn/\l+exp\(r— a)/c]}. Parame- 
ters a, c were taken the same as in |52j ]. namely a = 0.523 
fm and c chosen to satisfy i? rms = 0.836A 1 / 3 + 0.570 fm. 

Using the results of [H^ and this calculation we ob- 
tained in 0| the relative radiative corrections shown in 
Fig. El In order to include the interesting case Z = 55 
and to account for all values of 55 < Z < 90 we used 
interpolating formulae presented in [53, as well as data 
of 54] . The relative corrections for the Is and 2s lev- 
els are approximately the same size. This indicates that 
the radiative processes responsible for the correction take 
place at separations much smaller than the K-shell ra- 
dius, r <C (Zam)^ 1 , which is consistent with the as- 
sumption r < 1/m above. For these separations we can 
assume that, firstly, the screening of the nuclear Coulomb 
field in manyelectron atoms does not produce any sig- 
nificant effect, and, secondly, the relative corrections do 
not depend on the atomic energy level because for small 
separations all atomic nsi/2-wave functions exhibit sim- 
ilar behavior. These arguments remain valid for the pi i 2 
states as well, permitting us to presume that the results 
shown in Fig. 0] for the 2s!/2-levels and 2p 1 / 2 -levels of 
hydrogenlike ions remain valid for Siy 2 and P1/2 states of 
the valence electron in a manyelectron atom. We obtain 
the e-line radiative corrections for the PNC matrix ele- 
ment using the identity l|3.1[) that expresses them via the 
found corrections to the FNS energy shifts. The found 
PNC corrections, presented in Fig. 0]by the dotted line, 
are negative and large. For the 133 Cs atom the correction 
is -0.73(20)%, for Tl it is -1.6%. The result for 133 Cs 



is important because it indicated for the first time that 
the radiative corrections reconcile the experimental data 
of Wood et al [l]| with the standard model • 

A notable feature of the self-energy corrections pre- 
sented in Fig. 2] is their close similarity. The corrections 
to the FNS energy shifts for both si/2 and pi/ 2 states, 
as well as for the PNC amplitude, have large negative 
contributions from the QED self-energy, which monoton- 
ically and very smoothly increase with the nuclear charge. 
This observation brings us back to the point mentioned 
in section [J The large and negative contribution of the 
self-energy for heavy atoms is exactly what one should 
expect from it. Anything else would be a surprise. In 
order to verify this point one can try to apply the ap- 
proach outlined to some other problem, for example to 
the hyperfine interaction HFI which has been caref ully 
examined previously, see detailed discussions in 0, l47f 
and references therein [8l|. 

With this purpose one can try to derive the chiral in- 
variance identity for the HFI. The problem is that the 
HFI has a long-range tail ~ 1/r 3 that is prominent in 
the region where the chiral invariance is violated. How- 
ever, if one believes that convergence of the HFI matrix 
elements is fast enough, then the following equality 

^FNS, s ~ ^HFI, s ! (3.28) 

that was proposed in 0] should hold at least approx- 
imately. Here S' HF1 s is the radiative correction to the 
HFI for s 1 / 2 -levels, the primed notation indicates that 
the ^-independent Schwinger term a/(2n) should be ex- 
cluded. For heavy atoms this subtlety is not important, 
since the contribution from the Coulomb-induced correc- 
tions is much stronger than the Z-independent Schwinger 
term. Fig^ shows the e-line contribution to S' HFI s that 
was extracted in 0] from data of jH^ using interpola- 
tion for all values of Z considered there. It agrees semi- 
quantitatively with Ea. H3.28|) . the deviation is less than 
33 %. We observe again the same trend, the self-energy 
correction is large and negative, in accordance with the 
clear physical reasons mentioned in section [I] Another 
positive conclusion is that chiral invariance identities sim- 
ilar to H3.28|) can be applicable even for those potentials 
that spread out of the nuclear core. 

The approach outlined raises two points. Firstly, it is 
interesting to have an estimate for the self-energy correc- 
tion in simple terms, that would not appeal to sophisti- 
cated numerical calculations. Secondly, it is important 
to derive the self-energy corrections for lighter elements. 
Direct numerical calculations are difficult for them be- 
cause the error rapidly increases for lighter atoms. The 
main source for the error is the FNS energy shift for the 
Pi/2 state. The lightest element for which the self-energy 
corrections have been calculated is Cs, where the error 
for the FNS energy shift in the p\i 2 wave was ~ 100% 
[54l ] . Exactly this error results in the uncertainty of the 
above mentioned result —0.73(20)%. For heavier atoms 
the numerical errors of [52l | rapidly decrease. Therefore 
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the main error in our result for very heavy (Z > 80) 
atoms is associated with the validity of Ea. l|3.1l) . that 
was estimated above as a few percent. 

The estimate for the magnitude of the effect and cal- 
culations for lighter elements can be conveniently per- 
formed with the help of the Za-expansion considered in 
the next section. 



B. Perturbation theory in powers of Za 

Let us consider the perturbation theory for the e-line 
corrections. In the initial approximation one describes 
the electron using plane waves and consequently includ- 
ing in the perturbation theory three types of processes. 
Firstly, the PNC interaction, secondly, the e-line cor- 
rections, and, thirdly, the Coulomb interaction. Taking 
the nuclear Coulomb field as a perturbation one can for- 
mulate the perturbation theory in powers of the nuclear 
charge Z . The parameter that governs the correspond- 
ing expansion is a(Za) n , where n is the order of the 
perturbation theory. This statement follows from a con- 
ventional, well known result [6l| . In the Coulomb prob- 
lem the perturbation runs in powers of Ze 2 /hv, where v 
is a typical velocity of the electron. For short distances 
v ~ 1, thus Ze 2 /hv — Za. The type of expansion consid- 
ered typically includes some additional large logarithmic 
factors, as specified below. In the lowest order of per- 
turbation theory the simplest, linear in Z correction is 
oc Za 2 . Thus formulated problem requires calculation of 
the relevant coefficients that can be achieved using the 
corresponding Feynman diagrams. 

The Feynman diagrams for the e-line corrections in the 
first order of perturbation theory over the Coulomb field 
are presented in Fig. [5J The direct calculation of these 
diagrams for the PNC amplitude gives the following re- 
sult for the relative contribution of the e-line corrections 
to the PNC amplitude 

= -1.97 Za 2 . (3.29) 

It was first obtained in [5^ that had an amusing history 
|69|. The index 1 in the superscript of 1)3. 29[) reminds 
one that this is the linear in Z correction. For the 133 Cs 
atom Ea. (|3.29|) gives —0.6%, which is in line with the 
result —0.73(20)% that follows from the chiral invariance 
identity 13. 1|) as discussed in the previous section. An 
earlier attempt of Milstein and Sushkov gave the 
coefficient w +0.1, instead of -1.97 in Ea~ (|?H9"jl . This 
distinction had important implications. Based on their 
result the authors of claimed that there existed the 
contradiction between the experimental data of [TT| and 
the standard model, while Ea. (|3.29l) indicated that the 
experimental data agree with the model. Fortunately, 
the controversy was soon resolved by Milstein et al in 
the following Ref.j5|| that agreed with (|3.29|) and found 
also an analytical expression for the coefficient in this 



equation that reads -(23/4 - 4 log 2) = -1.970. This 
result has also been confirmed in recent Ref . |to| . 

At this point it is worth returning to the problem of 
gauge invariance. It is easy to demonstrate by conven- 
tional methods that the sum of all the diagrams in Fig. 
|3]is gauge invariant. This fact should be compared with 
the discussion of the gauge invariance in Section llH Al see 
after Ea. (|3.1|l . However, for practical applications one 
often fixes the gauge because this allows one to simplify 
lengthy analytical calculations. This is the way the calcu- 
lations were performed in [HH l56l Fol | . Since these works 
employed different gauges, the first one used the Feyn- 
man gauge while the second and third ones relied on the 
Fried- Yennie gauge (also called the Yennie gauge), their 
mutual agreement provides an additional helpful check of 
the validity of the final result Eq. (|3.29|l . 

Eg. 1)3.29(1 provides a simple transparent estimate for 
the contribution of the e-line corrections. One point to 
be noted is a strong dependence of the result on the nu- 
clear charge that makes the correction large for heavy 
atoms; remember the discussion of this point in Section 
[U This behavior of the self-energy correction differs qual- 
itatively from what we saw for the vacuum polarization in 
section^ where the higher-orders in the Za-expansion 
were found small. A notable feature in Ea. (|3.29|l is the 
large coefficient ~ —2.0 on the right-hand side. Naively 
one could expect this coefficient to be smaller, of the 
order of ~ 1/tt. It is interesting that a similar "numer- 
ical enhancement" occurs for the e-line radiative correc- 
tion for the energy shift that is due to the finite nuclear 
size (FNS). This correction was examined analytically in 
Refs.|67l l6Sl|. and more recently in Ref. [5(|. The result 
can be written as 

= -2.978 Za 2 , (3.30) 

where the analytical expression was derived in |56j |. We 
see that indeed, the e-line corrections to the FNS en- 
ergy shift are governed by the large coefficient ^ — 3.0 
in (GOP, similar to (|3~2^jl for the PNC amplitude. Fig. 
[5] examines this similarity in more detail. It shows data 
available for relative e-line corrections for the two prob- 
lems mentioned above, namely for the PNC amplitude 
and FNS energy shifts. The linear in Z approximations 
(|3.29|) . l|3.30|l (that are valid for sufficiently small values 
of Z) are compared in this figure with results that fol- 
low from numerical calculations of Ref. [52j . We observe 
a very close quantitative similarity between corrections 
to PNC and FNS. In both cases the linear approxima- 
tions H3.29fl and 13.30fl predict large negative corrections, 
which qualitatively agrees with results based on numeri- 
cal calculations for heavy atoms. The numerical validity 
of the linear approximations seems to be limited by the 
region below Z = 55, for higher Z they underestimate 
the effect. The Cs atom lies on the border, where results 
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of the small- iT and large- Z approaches agree reasonably 
well. 

Numerical data used in Eq. l|3.1|) incorporates an error 
that increases for smaller values of Z, see Tables III and 
IV of 52] . In order to reduce the impact of this error we 
can combine together Eqs. H3.1)l . (|3.29|l . One can approxi- 
mate all nonlinear terms omitted in (|3.29|) by the simplest 
quadratic function and choose the corresponding coeffi- 
cient to reproduce the results based on numerical data 
for very large Z, Z ~ 90, where the numerical errors are 
small. This approach gives the following interpolating 
formula [5{| for the corrections to the PNC amplitude 

^(e-nno,int) = _ L97 ^ ( X + 1>55 £ a ) _ ( 331 ) 

Fig. [fjjshows that the data available for large-Z and small- 
Z regions is very smooth. Therefore any reasonable in- 
terpolation would produce a pattern close to Ea. (|3.31|l . 
We can be certain therefore that l|3.31[l gives reliable 
numerical data. For the 133 Cs atom Eq. (|3.31[l predicts 
-0.9(1)%; compare this result with -0.6% of fHty and 
-0.73(20)% that follows from Eq.lJOJ H|. For the Tl 
atom Ea. (|3.31l) predicts —1.60% that agrees with the pre- 
diction of 1)3. 1[) discussed above (as it should, since the 
data of |52j] incorporates only very small errors for heavy 
atoms). 

It should be noted that Ea. (|3.1|l does not rely on the 
perturbation theory in powers of Za, effectively including 
all nonlinear corrections. This makes it accurate even at 
large values of the nuclear charge. It should therefore 
be considered formula that provides a convenient 
short-cut presentation of the nonlinear result. A similar 
formula was suggested in |5j| to reconcile the data for 
the e-line corrections to the energy shifts due to FNS for 
the Si/2 state. The linear in Z term of Ea. (|3.30|l and 
numerical data of |52j available for heavy atoms can be 
interpolated by 

^e-line, int) = _^ g7g ^2 ( j + g gg Za J (3 33) 

Fig. shows that this interpolation is reliable. This im- 
plies that for relatively light atoms Z ~ 55 — 75 numerical 
results for the Syii energy corrections should be amended 
by slightly larger values. 

An alternative method to include the nonlinear cor- 
rections was developed in Refs. [5rl l58| that relies on di- 
rect analytical calculations. This approach faces a diffi- 
culty because even in the next-to- leading order ~ Z 2 c? 
the calculations become complicated. Ref.^l showed, 
however, that in this order there appears a large log- 
arithmic factor that makes the correction proportional 
to - Z 2 o? lnl/mrN- Calculations with the logarith- 
mic accuracy, when some constant is neglected, are much 
more feasible. The elegant result for the thus calculated 
second-order correction (index 2 in the superscript) was 
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where b = exp[l/(2 7 ) - C- 5/6]. The fina l resu lt of ^ 
is given by the sum of the linear term Eq. (|3.29|l and the 
second order correction Ea. (|3.33|l . It is compared in Fig. 
[3 with predictions of Eqs. (|3.1|l and l|3.31|l . For the sake 
of completeness the figure also shows the previous results 
of Ref . • For all values of the nuclear charges the re- 
sults of |56j are close to our results derived from Eas. (|3.1|) 
and (|3.31|) . as was emphasized in |^. Recently Milstein 
et al |58j refined their arguments, demonstrating that the 
logarithmic factor similar to the one in Ea. H3.33[) exists 
in all higher order terms of the Za-expansion. Ref. [5^ 
found the corresponding contribution of the higher-order 
terms numerically. This latest calculation, also shown in 
FigQ is even closer to our results given by Eas. (|3.1|l and 
(|3.31(l . The deviation from (|3.31l) is below 9% for a wide 
range of atomic charges. Ref. £3 claims their error to 
be of the order of ~ 5%; our estimate for the accuracy of 
Ea. (|3.1|) . which limits the error at large Z, is a few per- 
cent, see Section llll Al Within these errors the results of 
the two groups completely agree. The fact that the re- 
maining discrepancy increases for heavy atoms indicates 
probably that it is due to terms of the Za-expansion still 
unaccounted for in (5^. Overall, Fig|7|shows good agree- 
ment of the recent results obtained by the two groups 
using different approaches. The sharp contradiction that 
existed during the initial stages of this research (com- 
pare the dashed-dotted line with the thick-dotted line) 
makes the latest convergence even more satisfying and 
trustworthy . 

There is an interesting physical link between the chi- 
ral invariance identity of Ea. (|3.1() and the perturbation 
theory approach of Eas. (|3.29|l . (|3.33|) . Remember that in 
deriving Ea. (|2.3l) in the previous section we focused our 
attention at one stage of the consideration at small dis- 
tances r < Zot/m. The claim was that if this region is 
proved to comply with Ea. (|2.3|l . then the outer region 
r > Za/m would inevitably comply with it as well, see 
discussion after Ea. l3.24|) . A similar philosophy lies be- 
hind the logarithmic approximation in Eq. (|3.33|) . The 
integration that leads to the logarithmic function in this 
equation is saturated at small distances. The omitted 
constant is related to larger distances, but it is of lesser 
importance and can, in the simplest approximation, be 
neglected. Thus in both approaches the region of small 
distances proves to be the most important. 

This link between the chiral invariance identity and the 
perturbation theory shows that Ea. l|3.1|l should be valid 
when one restricts consideration to the logarithmic ap- 
proximation. This means that the term oc Z 2 o? InmrN 
should exist not only in the PNC problem, see Eq. l|3.33|l . 
but in the FNS problem for the Sj/2 and P1/2 partial 
waves as well. Ea. l|3.1ll predicts a linear relation be- 
tween the three coefficients of these three terms. Ref. 15a 
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demonstrated that the three coefficients mentioned are 
all equal, in compliance with Eci. (|3.1|) . This fact can be 
considered either as a consequence of Ea. (|3.1|) . or, al- 
ternatively, as its independent verification, being fruitful 
either way. 

The picture of reliability presented above based on the 
complete accord of all available data for the self-energy 
corrections in atoms has recently been put to test again 
by Milstein et al. They have claimed in Ref. j5(| that for 
small values of the nuclear charge Z ~ 1 the self-energy 
corrections for the energy shift induced by the FNS in the 
P1/2 partial wave are positive and very large, being two 
orders of magnitude greater than anticipated previously 
for the hydrogen atom. For the p 3 / 2 state the proposed 
enhancement is even stronger, four orders of magnitude. 
This conclusion, if correct, could be related to the PNC 
problem for light elements through the equality of 
mentioned above that binds together the radiative correc- 
tions for the PNC and FNS problems. It needs, however, 
to be p ointed out that a questionable approach is applied 
in 59] . It uses the small- momenta expansion for the ver- 
tex operator. When this technique is used for the short- 
range potential there arises a contradiction because the 
short-range potential needs large momenta. Attempts to 
use a similar line of arguments for calculations of radia- 
tive corrections have been made previously, but it was 
recognized that they strongly overestimate the effect, by 
orders of magnitude. A necessity to be very careful, to 
avoid using the small-momenta asymptotic dealing with 
the short-range potentials, was clearly stated long time 
ago by Lepage et al [6(J • One is forced to presume there- 
fore that conclusions of |59j| are not convincing (probably 
they should be even called preliminary, if one uses the 
latter term as in Ref. 01)- Fortunately, these results do 
not influence either the above mentioned good numeri- 
cal agreement that was achieved for heavy atoms, or the 
analytical agreement between the two groups within the 
Za-expansion that defines the PNC in light elements. 

Discussions in this section are almost entirely devoted 
to analytical methods, the only exception was made for 
the numerical data available for corrections to the energy 
shifts due to the FNS. The main reason for this is that 
direct numerical calculations for the self-energy correc- 
tions for the PNC amplitude have not yet been reported. 
Presumably the difficulties in this approach are signifi- 
cant, though not being experts in this area we would not 
speculate on the details. We wish to mention, however, 
that developments in this direction are quite desirable 
due to several reasons. One of them, purely theoretical, 
has come into existence only recently. Numerical studies 
can determine independently the accuracy of the chiral 
identity (|3.1[1 . Knowing the limitations of this relation in 
the PNC problem, one would be able to estimate the ac- 
curacy of other relations of this type that can be derived 
and used in similar problems in the future. 

Summarizing, we discussed above two different meth- 
ods for calculation of the self-energy corrections, one re- 
lated to the chiral invariance identity (|3.1[) . the other one 



based on the Za-expansion. They give close results prov- 
ing that the self-energy corrections to the PNC amplitude 
are large and negative. Most importantly, this fact brings 
the experimental data for the 6s — 7s PNC amplitude in 
133 Cs into agreement with the standard model. 

IV. COMPARISON WITH EXPERIMENTAL 
DATA 

Fig. summarizes the data for the QED radiative cor- 
rections presenting the vacuum polarization correction 
4^0 calculated with the help of Ea. $FFty - $FZty and 
the self-energy (plus vertex) radiative correction <5pNC 
described by Eq. I|3.31|) , as well as the total radiative cor- 
rection 

,(tot) _ .(VP) .(e-line) -v 

°pnc — °pnc ' °pnc • v^- 1 ; 

We see that both the vacuum polarization and the 
self-energy corrections rise quickly with the nuclear 
charge, but their opposite signs make the total correc- 
tion smoother. We find <5pNC = —0.54% for the Cs atom 
and 4nc = -0-70% for Pb,Tl and Bi. 

Let us discuss the implications of these results, first 
for the 6s - 7s PNC amplitude in 133 Cs. The standard 
model value for the nuclear weak charge for Cs |7l| is 

QwC 33 Cs) = -73.09 ± (0.03) . (4.2) 

Ref. [37j refined previous calculations of Ref. ex- 
tracting from the experimental PNC amplitude of Ref. 
[Tl| the weak charge 

Q (C +B+ N) (13 3 Cs)= (43) 

-72.45 ±(0.29) cxpt ±(0.36)theor, 

with the theoretical error 0.5%. This value includes the 
correlation and the Breit corrections, as well as the neu- 
tron skin corrections (superscript C+B+N) , but does not 
take into account the radiative corrections 83] . The neu- 
tron nuclear radius is slightly larger than its proton ra- 
dius. This provides a possibility for the electron to in- 
teract with neutrons in the region outside of the proton 
core, where the wave function is slightly smaller than 
inside the nucleus. The corresponding reduction of the 
PNC amplitude is called the neutron skin correction in 
PNC. According to estimates of Derevianko the neu- 
tron skin correction is approximately —0.2% for Cs. The 
calculations of Ref. 37J took this effect into account. 

Eq.lO is consistent with Q w ( 133 Gs) = -72.21 ± 
(0.28) OX pt ± (0.34) t hoor that was adopted in fill by tak- 
ing the average of the results of Refs. |3%l3al35| . and ac- 
cepting the theoretical error 0.4% proposed in |34|. This 
value for the weak charge includes the vacuum polariza- 
tion correction 0.4%. Extracting it and taking into ac- 
count —0.2% for the nuclear skin correction one obtains 
-72.35 that is close to Eq.O- 



17 



The radiative corrections derived from results pre- 
sented in Fig. [8] are — 0.54 ± (0.10)%; the error reflects 
the uncertainty of the self-energy radiative correction. 
Eci. (|4.3|) combined with the QED radiative corrections 
gives 

Q Total ( 133 Cs) = (4 4) 

-72.84 ±(0.29) expt ±(0.36) theor . 

where Total = Correlations + Breit corrections + Neu- 
tron skin correction + QED radiative corrections. The 
agreement with the standard model given in Ea. (|4.2|l is 
good. It is so good, in fact, that we hasten to remind the 
reader a historical aspect of the problem. The theory has 
come up with Eq. 14. 4|) after a turbulent period of research 
during which the experimental data were widely antici- 
pated to be in contradiction with the standard model. 
This makes the found "unexpected" agreement more ob- 
jective. 

Let us consider now the case of the thallium atom that 
was studied experimentally in H3. Il3l Ill HM llH Il7 | . the 

calculations were performed in |74l that give close 
results for the many-electron correlation. The corrections 
to the result of should include contributions of the 
Breit interaction -0.98% jz2]. As was mentioned, the 
QED radiative corrections for Tl give —0.70%. If we 
adopt the neutron skin correction —0.2%, then the total 
theoretical result reads 

R = Im(£pNc/Ml) = -15.64(-Qi V /A r )T0~ 8 . (4.5) 

Comparing with the measured value of Ref.^i| that pre- 
dicts R = —14.68(17) • 10~ 8 , we obtain the weak charge 
for 205 Tl 

QwC° 5 Tl) = -H6.4 ± (1.3) cxpt ± (3.4) thoor (4.6) 
in good agreement with the standard model prediction 

Qw( 205 T\) = -116.7(1) . (4.7) 

V. CONCLUSION 

The QED corrections for parity nonconservation have 
recently emerged as an important ingredient in the anal- 
yses of modern experimental data. The improvement in 
the experimental accuracy, particularly for the 6s — 7s 
PNC amplitude in 133 Cs, and the progress of atomic 



structure calculations revealed a possible gap between 
the experiment data and atomic structure calculations 
on one side and the standard model on the other. 

The QED corrections, which embrace the Breit cor- 
rections and the radiative corrections, reconcile the ex- 
perimental data with the standard model. The area has 
evolved and developed very rapidly. Over the past couple 
of years, even over a year, the situation has changed dra- 
matically. A short time ago one could be seriously con- 
templating the possibility to modify the standard model 
to accommodate a possible deviation of atomic experi- 
mental data with the standard model. Presently it has 
become clear that modifications of the model are to be 
postponed, one simply needs to calculate everything ac- 
curately. Everything here is essential. Different parts of 
the QED corrections are all very large and have different 
signs. One must be certain that all parts are included 
and accounted for properly. That is why we spent some 
time and space above presenting and comparing all avail- 
able results for the radiative corrections. The conclusion 
is that all data for the radiative corrections fit together 
well. For each possible correction several different meth- 
ods of calculation were employed by different groups with 
the same final conclusions. The self-energy corrections, 
that proved to be by far the most difficult ones, are known 
presently with an uncertainty of better than 10 %, which 
is sufficient for the present day experimental accuracy. 
One can contemplate a significant reduction of this er- 
ror, if required. 

The calculations of the radiative corrections for PNC 
were based on vast expertise accumulated in other related 
research areas that include the Lamb shift, the hyperfine 
interaction, the nuclear finite-size correction. In turn, 
some methods of calculations that have recently evolved 
in relation to PNC may find applications in other areas. 
One of them is the chiral invariance identity Eq. l|3.1fl . 
that was originally applied to express radiative correc- 
tions to PNC in terms of the corrections to energy shifts 
due to the finite nuclear size. Similar identities can be 
derived for other related problems. One can also mention 
a simple convenient expression that gives the variation of 
the electron wave function at the atomic nucleus in terms 
of the two first moments of the perturbative potential 
Eci.(l2~31). 

In conclusion, large QED radiative corrections to the 
parity nonconservation amplitude in heavy atoms rec- 
oncile the accurate atomic experimental data on parity 
nonconservation for the 6s — 7s PNC amplitude in 133 Cs 
of Wood et al 0] with the standard model. 

Discussions with D.Budker, J.S.M.Ginges, W.Greiner, 
G.F.Gribakin, M.G.Kozlov, V.M.Shabaev, and 
A.V.Solov'yov are greatly appreciated. This work 
was supported by the Australian Research Council. 
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FIG. 1: Relative contribution (%) of the QED vacuum po- 
larization to the PNC amplitude versus the nuclear charge. 
Thick line - total correction predicted by Eq. l2.21ll . dashed 
line - prediction of Eq. (12.2211 that takes into account only 
events outside the nucleus, dotted line - results based on the 
expansion Eg. 12. 251 . 



FIG. 4: The relative radiative corrections (%) induced by 
the e-line corrections (self-energy plus vertex) of Fig|5] Thick, 
thin, and long-dashed lines - corrections for FNS energy shifts 
for ls-i/2, and 2p 1 / 2 levels respectively extracted from 

[5"H|: short-dashed line - prediction of Eg. 13. Ill for the PNC 
matrix element calculated in |48| , dashed-dotted line - correc- 
tion to the hyperfine interaction |4ql . 



FIG. 2: The Feynman diagrams for QED vertex (a) and 
self-energy (b),(c) corrections, called e-line corrections in the 
text. They are used in the text to describe either the PNC 
amplitude or the energy shifts induced by the finite nuclear 
size. The thick dot represents in these two cases either the 
vertex induced by the PNC Hamiltonian, or the finite nuclear 
size (FNS) potential that takes into account the spread of the 
nuclear charge inside the nucleus. The thick line describes the 
electron propagation in the atomic field. 



FIG. 5: The lowest order Feynman diagrams for the QED 
self-energy and vertex corrections, called e-line corrections in 
the text, to the PNC matrix element. The thin line describes 
propagation of a free electron. For each diagram one of the 
wavy legs shows the Coulomb interaction with the nucleus, 
another one - the weak PNC interaction with the nucleus 
(which was shown by a large dot in Fig. |5J). Each diagram 
represents all possible Feynman diagrams with the given topo- 
logical structure. 



FIG. 3: The S1/2 and pi/2 wave functions for the valence 
electron in 133 Cs (arbitrary units) versus the distance r/rc, 
where rc = h/mc is the Compton radius. The functions are 
normalized in such a way that the const in both sets of bound- 
ary conditions in Eg. 13. 1611 takes the same value, const — 1. 
Solid and dotted lines - / s ,i/2 and <? Pi i/2 wave functions that 
are large for this normalization, thick dashed and thin dashed 
lines —g s> i/2 and / p ,i/2 wave functions that are much smaller, 
(a) the region of large distances r ~ rc shows good coinci- 
dence of / s ,i/2 with g p ,i/2, (b) the region in the vicinity of the 
nucleus, the nuclear radius is rjj — 0.016 rc — 6.1 frn, where 
— 9s, 1/2 coincides with / Pl i/2- The shown results illustrate 
Eo. l3.6H based on chiral invariance. 



FIG. 6: Thick dotted, thin dotted and thick solid lines show 
e-line corrections (self-energy plus vertex) to the PNC ampli- 
tude predicted by Eq . (12.311 of 48], by the linear approxima- 
tion Eo. 13.29H of |55|. and by the interpolating Eo. 13.3H of 
|55^l (all same as in Fig. 0. The thin dashed line, thick dashed 
line, and the dashed-dotted line - the e-line corrections to the 
FNS ener gy shifts in the linear approximation 13.3U1 . from 
results of |5lll . and the interpolating formula 13.32H . Note the 
similarity between the PNC and FNS problems. 
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FIG. 7: Relative contribution of the e-line corrections (self- 
energy plus vertex) to the PNC amplitude ( % ) versus the 
nuclear charge. Thick dotted, thin dotted and thick solid lines 
show predictions of Eq. 12.311 of |4g|. of the linear approxi- 
mation i3~%?t of [13, and the interpolating Eq. ll33Tl of j5^| 
respectively, as in Fig. [(J Dashed-dotted line, thin dashed line 
and thick dashed line - results of Milstein et al from |43 . l56ll5S| 
respectively. Note the convergence of the latest results of the 
two groups, which strongly deviated at the initial stage. 



FIG. 8: QED radiative corrections to the PNC amplitude (%) 
versus the nuclear charge. The dotted line - vacuum polar- 
ization Eq. 12.21L thin line - self energy correction Eq. 13.31L 
thick line - total correction Eq. I4.1L 
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